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Abstract 

Using the framework for solving the spherically symmetric field 
equations in observational coordinates given in Araujo & Stoeger 
(1999), their formulation and solution in the perturbed FLRW spheri- 
cally symmetric case with observational data representing galaxy red- 
shifts, number counts and observer area distances, both as functions 
of redshift on our past light cone, are presented. The importance of 
the central conditions, those which must hold on our world line C, 
is emphasized. In detailing the solution for these perturbations, we 
discuss the gauge problem and its resolution in this context, as well 
as how errors and gaps in the data are propagated together with the 



genuine perturbations. This will provide guidance for solving, and in- 
terpreting the solutions of the more complicated general perturbation 
problem with observational data on our past light cone. 



1 Introduction 



In two recent papers ( |[Araujo and Stoeger (1999)1 , |Araujo et al. ( 2001 )|| ) we 



demonstrated in detail how to solve exactly the Einstein field equations for 
dust in observational coordinates with cosmological data function represent- 
ing galaxy redshifts, and observer area distances and galaxy number counts 
as functions of redshift. These data are given, not on a space-like surface of 
constant time, but rather on our past light cone C~(po), which is centered 
at our observational position po "here and now" on our world line C. These 
results demonstrate how cosmologically relevant astronomical data can be 
used to determine the space-time structure of the universe - the cosmolog- 
ical model which best fits our universe. This has been the aim of a series 
of papers going back to the Physics Reports paper by Ellis et al, (1985) 
The motivation and history of this "observational cosmology program" is 
summarized in Araujo & Stoeger (1999). 

The primary aim of this program is to strengthen the connections be- 
tween astronomical observations and cosmological theory. We do this by 
allowing observational data to determine the geometry of spacetime as much 
as possible, without relying on a priori assumptions more than is necessary 
or justified. Basically, we want to find out not only how far our observ- 
able universe is from being isotropic and spatially homogeneous (that, is 
describable by a Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmolog- 
ical model) on various length scales, but also to give a dynamic account of 
those deviations ( ||Stoeger et al. (1992a) 



Although there are strong indications, especially from the character of the 
cosmic microwave background radiation (CMWBR), ( fStoeger et al. (1995 



[Stoeger et al. (1997)|| ) that our universe is very close to being FLRW on the 
largest length scales, it is very clear that, since it is lumpy on small and 
intermediate length scales, it is not exactly FLRW. In light of this it is very 
important to determine, from the data we in principle have available, how 
the universe deviates from FLRW on the largest scales, and how those de- 
viations grow or damp as we move off our past light cone into the past, or 
into the future, as well as the degree to which the errors and gaps in the 
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data induce imprecisions in the cosmological model we adopt. Genuine per- 
turbations from FLRW and errors in the data will be propagated together 
( [[Stoeger et al. (1992a)|| ). In practice, they will have to be separated out 
on C~(po) and then followed individually in time. Both of these goals are 
pursued in this paper and in a subsequent paper, where we show in detail 
how to solve the perturbed field equations with cosmological data on C~(p ) 
with a best-fit FLRW model as the background, or zeroth- order, solution. In 
this paper we limit ourselves to spherically symmetric perturbations, using 
the framework we developed in our recent exact spherically symmetric treat- 
ment ( [|Araujo and Stoeger (1999~| ). In the subsequent paper we shall tackle 
the case of general perturbations. These two cases were previously studied 
in two much earlier papers, ([ [Stoeger et al. (1992a)| , |Stoeger et al. (1992b)| |) 
but unsuccessfully, due to an important error which led to the severe over- 
restriction of the solutions (see Araiijo & Stoeger (1999) and references 
therein) . 

What is the relevance of constructing spherically symmetric perturbed 
solutions to FLRW, when we could use the same data to find the exact spheri- 
cally symmetric solution ( [[Araujo and Stoeger (1999)|| )? It is obviously very 
important to understand observationally based perturbations to FLRW in 
their own right. Having a secure approach for determining exact spherically 
symmetric solutions enables us to study and interpret spherically symmetric 
perturbations very simply and carefully, including the gauge problem (see 
Section 4) and the propagation of errors and gaps in the data, as mentioned 
above, before going on to solve and and interpret the more complicated gen- 
eral perturbation case in this observational context. Finally, studying obser- 
vationally determined spherically symmetric perturbations will help us define 
criteria for the use of perturbation theory itself - or, equivalently, for defin- 
ing a class of 'almost FLRW models and for determining whether or not our 
observable universe is 'almost FLRW,' and on what length scale, in a mean- 
ingful and rigorous way. Certainly, it is not 'almost FLRW on small and 
many intermediate length scales, but there is strong evidence from CMWBR 
measurements that it is on the largest scales, as we indicated above. At 
what length scale can we begin to use 'almost FLRW models to describe the 
universe? Using this approach, we should be able eventually to answer this 
important question. Further discussion of the philosophy of this approach is 
given in Stoeger et al, (1992a). 

In Section 2, we briefly summarize our characterization of observational 
coordinates, as well as the metric and observational relations in observational 
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coordinates, within the framework of spherically symmetric deviations from 
FLRW data. In Section 3, we give the exact spherically symmetric field 
equations in observational coordinates, and in section 4, we discuss gauge 
problem as it relates to our observationally based cosmological perturbations 
from FLRW. We describe in Section 5 the first step of the integration proce- 
dure, which is the solution of the perturbed null Raychaudhuri equation to 
determine the redshift function of the null radial coordinate y, that 

is to find z = z(y). Finally, in Section 6, we complete the integration proce- 
dure by explicitly solving for all the spherically symmetric perturbations in 
terms of the data functions of z\{y), the non-FLRW component of z(y), and 
we briefly discuss the application of these results to the issues we have just 
highlighted. 

2 Coordinates, Metric and Observational 
Relations 

We are using observational coordinates (which were first suggested by Temple 
(1938)). As described by Ellis et al, (1985) these coordinates x % = {w, y, 9, 0} 
are centered on the observer's world line C and defined in the following way: 

(i) w is constant on each past light cone along C, with u a d a w > along 
C, where u a is the 4- velocity of matter (u a u a = — 1). In other words, each 
w = constant specifies a past light cone along C. Our past light cone is 
designated as w = wo- 

(ii) y is the null radial coordinate. It measures distance down the null 
geodesies - with affine parameter v - generating each past light cone centered 
on C. y = on C and dy/du > on each null cone - so that y increases as 
one moves down a past light cone away from C. 

(iii) 9 and <fi are the latitude and longitude of observation, respectively - 
spherical coordinates based on a parallelly propagated orthonormal tetrad 
along C, and defined away from C by k a d a 9 = k a d a (j) = 0, where k a is the 
past-directed wave vector of photons {k a k a = 0). 

There are certain freedoms in the specification of these observational coor- 
dinates. In w there is the remaining freedom to specify w along our world 
line C. Once specified there it is fixed for all other world lines. There is 
considerable freedom in the choice of y - there are a large variety of possible 
choices for this coordinate - the affine parameter, the redshift z, the observer 
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area distance C(w,y) itself, which we oftentimes write as r (z) or r (y) on 
our past light cone w = Wq. We normally choose y to be comoving with the 
fluid, that is u a d a y = 0. Once we have made this choice, there is still a little 
bit of freedom left in y, which we shall use below. The freedom in the 9 and 
(j) coordinates corresponds to a rigid rotation of orthonormal tetrad at one 
point, say p , on our world line C. They are just spherical coordinates on the 
celestial sphere with respect to the (physically non-rotating) reference frame 
of the orthonormal tetrad. 

Since we are using the best-fit FLRW model as our background space- 
time the most general pertubed metric in observational coordinates takes the 
form: 



V 



-R 2 + Z 2 


R 2 + f3 2 


v 2 


V3 


R 2 + (3 2 











V 2 





R 2 f 2 + h 22 


h 23 


v 3 





h 23 


R 2 f 2 sin 2 9 + h 33 



where R 2 , R 2 f 2 and R 2 f 2 sin 2 9 are the FLRW values of the respective metric 
components, taken here as zeroth-order terms, and all the other terms are 
the non-FLRW perturbations in the sense described in the introduction. / 
is given by 



sin y k = 1 (closed) 
y k = (flat) 

sinhy k = — 1 (open). 



(2) 



In this paper we are concerned with spherical perturbations. The metric 
then is a lot simplified with 



v 2 = v 3 = 



(3) 



h 



23 







(4) 



^22 = /i33(sin6>) 2 . (5) 

Therefore, the spherically symmetric perturbation problem is reduced to 
solving the linearized field equations for j3 2 (w,y), Z 2 (w,y) and h 22 (w,y). 
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The remaining coordinate freedom which preserves the observational form 
of the metric is a scaling of w and of y: 



w _> w = w(w) , y^y = y(y) ^ ^ . (6) 

The first, as we mentioned above, corresponds to a freedom to choose w 
as any time parameter we wish along C, along our world line at y — 0. This 
is effected in this case by choosing (— goo(u>, 0)) 1//2 . The second corresponds 
to the freedom to choose y as any null distance parameter on an initial light 
cone - typically our light cone at w = wq. Then that choice is effectively 
dragged onto other light cones by the fluid flow - y is comoving with the 
fluid 4-velocity, as we have already indicated. We shall use this freedom to 
choose y by setting: 

(3 2 (w Q ,y) = -Z 2 (w ,y) (7) 

We note that the first attempt to solve this problem by [Stoeger et al. 
(1992a)] is in error as it relies on the assumption that (3 2 {w,y) = —Z 2 (w,y), 
which is too restrictive when y is comoving. 

Before we proceed we write down our observational quantities in terms 
of the coordinate metric components given above: 

(i) Redshift. The redshift z at time Wq on C for a comoving source a null 
radial distance y down C~(po) is given by 

l + z = A [R 2 (w ,y) - Z 2 (w ,y)]- l/2 , (8) 

where, A = [R 2 (wo,0) — Z 2 (wo, 0)] 1 / 2 . Expanding (||) to first order using 
the binomial theorem yields 



1 + 2 2 R(w a ,y) 

This is just the observed redshift, which is directly determined by source 
spectra, once they are corrected for the Doppler shift due to local motions. 

(ii) Observer Area Distance. The observer area distance, often written 
as r , measured at time il>o on C for a source at a null radial distance y is 
given by the expression r 4 sin 2 8 = det(gjj), I, J ranging over the values 2,3 
([ |Ellis ct al. (198"5]| ). In our case, it takes the form: 
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ro = [R 2 (w ,y)f 2 +h 22 (w ,y)} 1 / 2 = R(w ,y)f(y)+ ^ 22 +.... (10) 

2R(w ,y)f{y) 

(iii) Galaxy Number Counts. The number of galaxies counted by a cen- 
tral observer out to a null radial distance y is given by (see [Ellis ( 1971 )| , 
Ellis et al. (19851 ): 

N(y) = 4vr f\{w^y) m - 1 (R 2 (y) - Z 2 (y)) 1/2 r (y) 2 dy (11) 

where /i is the mass-energy density, m is the average galaxy mass, R = 
R(w Q ,y) and Zq = Z 2 (w ,y). Then the total energy density can be written 

as 

dz 1 

fi{w , y)=m n(w , y) = M Q {z) — — ^ _ ^ /2 (12) 

where n(wo,y) is the number density of sources at (wo,y), and where 

„ , / s ml 1 dN . „ s 

Here dQ is the solid angle over which sources are counted, and J is the 
completeness of the galaxy count, that is, the fraction of sources in the 
volume that are counted is J. The effects of dark matter in biasing the 
galactic distribution may be incorporated via J. In particular, strong bi- 
asing is needed if the number counts have a fractal behaviour on local scales 
( [Humphreys et al. (19981 ). 



3 The spherically symmetric field equations 
in observational coordinates 

In this section we present the exact spherically symmetric field equations 
([ |Araujo and Stoeger (1999~| ) in observational coordinates for the case of 



dust. The perturbed spherically symmetric field equations can be constructed 
from them, as we show in the sequel. 

In observational coordinates the spherically symmetric metric takes the 
general form: 

ds 2 = -A(w, yfdw 2 + 2A(w, y)B(w, y)dwdy + C{w, y) 2 dtt 2 . (14) 
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The central conditions for the metric variables A(w,y), B(w,y) and 
C(w,y) in ([14]) - that is, their proper behavior as they approach y = 
are: 



as y : A(w, y) -> A(w, 0) ^ 0, 
B(w,y)^ %,0)^0, 
C(tu,y)-B(w,0)2/ = 0, (15) 
C„(u;,jO->B(u;,0). 

As pointed out in Araujo & Stoeger (1999), in this exact case the freedom 
to choose the coordinate y is used by setting 

A{w ,y)=B{w ,y) (16) 
Also the redshift at time wo is given by 

! + ^ = ^A^ ( 17 ) 
A(w ,y) 

In our present notation A, B and C are identified as 

A = [R 2 (w,y)-Z 2 (w,y)] 1 / 2 (18) 

= [fl>,y) + /?>,y)] f , 

" -Z> )2 /)]V2 1 1 

C = [R 2 (w,y)f 2 (y) + h 22 (w,y)] 1 / 2 . (20) 

As firstly pointed out by Stoeger et ai, (1992c) the general spherically 
symmetric field equations for dust in terms of the metric ( |14|) split in two 
sets of radial and time equations. 

The time independent equations are 



CJ C A 1 /^^x 

? = c~ A ( ' 

B B A , r, 1 . . . 

— = 2uA 2 - -uA 2 , 22 

B BA 2 P ' V ; 

where an overdot indicates 9/<9u> and the quantities uj and /x above are given 
by 
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H(w,y) = /x„(y) B-\w,y) C- 2 (w,y) (23) 

/ \ My) i c c ifC\ 2 

u{w/y) = CH^ = -^ + ACBC + 2{Bc) ' (24) 

where a prime indicates d/dy, // again is the relativistic mass-energy density 

of the dust and u is a quantity closely related to /x (see equation Q2DD 
below). Both u; and /x are specified by data on our past light cone. 
The radial equations are 



C" 

~c 




(25) 



<{y) = -^o(y) + ^ ) (26) 

The contracted Bianchi identities yields 

A' + 5 = 0. (28) 
From EBJ we see that there is a naturally defined "potential" 



W(y) = j + (29) 

that simplifies the integration procedure since it implies that the time- 
derivative equations (PH and |22|) are not necessary to obtain the solution 
becoming consistency conditions ( [|Araujo and Stocgcr (1999)|| ). 

Stoeger et ai, (1992c) and Maartens et ai, (1996) have shown that equa- 



tions (|24|) and (|29|) can be transformed into equations for A and B, thus 



reducing the problem to determining C: 



A 
B 



C 

[W 2 - 1 -2cj /C] 1 /2 

C_ 

W- [W 2 - 1-2cj /C] 1 /2- 



(30) 
(31) 
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Now, as pointed out by Araujo and Stoeger (1999), if we temporarily 
choose y = z, which is always legitimate, knowing C(w ,z) from the data, 
lvo(z) can be obtained from equation (|3T| ) and has the form 



CW 5 



U 



1 



1 



WA 



with y 



z, and we have used ([16]) and ( 17|) to write 
B(w ,y)=A(w ,0)/(l + z). 



(32) 



(33) 



4 The Gauge Problem 

Before we proceed to integrate the perturbed spherically symmetric field 
equations we address in this section the gauge problem, a problem that arises 
when one considers perturbed models in General Relativity due to the fact 
that there is no invariant way to define a background space-time given the 
actual space-time. 

Perturbations are usually discussed by assuming that the space-time (M, 
g) admits a family of coordinate systems in which the metric tensor compo- 
nents can be written as 



9„u(x a , v) =V M + Vh^{x a ) + -^k^) + yV%»(x a ) + ■■■ (34) 

where g, h, k,l . . ., are of the same magnitude and rj <C 1. Here we are 
comparing a point in the space-time (M, g) with coordinates x a to a point 

o 

in another space-time (background) (M ,g) with the same coordinates. 

Following Stewart and Walker (1974), we adopt in the present discussion 
the point of view that the perturbed space-time (M, g) is to be thought of as 

o 

the result of some slight changes made in the background space-time (M , 9). 
This approach emphasizes the notion that the background and the perturbed 
space-times are considered as distinct objects. One also requires the conti- 

o 

nuity of the perturbation by assuming that (M , 9) and (M, g) are connected 
by a path in the space of space-times. Therefore, since the perturbation is 
small, we consider sequences of diffeomorphic space-times depending on the 
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continuously varying parameter r\. As we shall see below, the concepts of co- 
ordinate and gauge transformations, as well as the linearization of the field 
equations, are made more transparent when viewed in this light. 

In order to formalise these ideas (IIGeroch (1969)11), we assume there ex- 



ists a smooth, Hausdorff, five dimensional manifold M in which the space- 
times under consideration are smooth, nonintersecting, properly-embedded 
four- dimensional submanifolds. In other words, we are considering a one- 
parameter family of space-times [M v ,g(rj)] embedded in five-dimensional 
manifold M. The parametrisation is chosen in such a way that the back- 
ground space-time is given by r\ = 0. 

We now consider any smooth, nowhere-vanishing vector field v on M 
which is transversal (nowhere tangent) to the M n . This vector field naturally 
defines a map, : M — > M v , between the background and the perturbed 
space-times which identifies points lying on the same integral curve of v. This 
map is called the identification map because it actually defines when points 
in different space-times are to be regarded as the same. A choice of vector 
field v is called a choice of identification gauge. 

Having made these remarks, we now recognise (0) as arising from a 
choice of identification gauge such that : M — > M v maps points having 
the same coordinates. 

We consider two kinds of coordinate transformations within the family of 
coordinate systems in which ([34]) holds. Under the first we make the same 
coordinate changes on both M and M v and keep identifying points with 



the same coordinates, i.e. all terms in (34) transform as tensors. Of course, 
h^i k^ u , ■ ■ ■ are not tensors by themselves but rather parts of g^ u . However, in 
analogy to the weak field approximation in which the background space-time 
is flat, the idea is to regard h^, k^, ... as tensor fields on the background 
space-time. This transformation is therefore, the strong field analogue to 
the background Lorentz transformations of the weak field approximation. 
Clearly, it does not affect the ordering scheme of ( p4] ) and it is obvious that 
we are still in the same identification gauge. So, the coordinate freedoms 
mentioned in section 2 fall into this class. 

Under the second kind of coordinate transformations we make different 
coordinate changes on Mo and M v and still identify points with the same 
coordinates. To illustrate this important point, suppose we make a small 
coordinate transformation on Mq of the form 

x » _> = x^ + r]C (35) 
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where ^ is an arbitrary vector field on M . On M eta we do not make 
any coordinate transformation. Then, neglecting terms of 0(r] 2 ), we find 
( ||Landau and Lifshits (19751 ) 

hpv - h^ v = 2£ (At;i ,) = g)^ (36) 

where £g is the Lie derivative with respect to the vector field £. We observe 
from this example that the whole effect of the coordinate transformation is ex- 
pressed in h^. This coordinate transformation also does not alter the order- 
ing scheme given by ([34]), so we are still in the assumed family of coordinate 
systems. However, it shows that h^ v is not fixed uniquely, or equivalently, 
that the background space-time is only fixed up to small transformations of 
this kind. These transformations are called gauge transformations. Under a 
gauge transformation all quantities dependent upon h^ v in general also un- 
dergo transformations. It is clear from ([36]) that if a quantity vanishes in the 
background space-time, then it is gauge invariant to all orders. 

In our case here we are using a best-fit FLRW model as our background 
space-time, and, as we have pointed out above, the further coordinate trans- 
formations allowed, respecting this choice, are those we have labelled as type 
1. As we shall see below (Section 6), given the data, and the coordinate 
choices we have insisted upon, along with the central conditions, equations 



( 15|) , there is practically no coordinate freedom left. Essentially, our obser- 
vational frame of reference and the coordinates we use to express that along 
with the choice of the FLRW background by best-fitting methods is equiva- 
lent to choosing the gauge. From this point of view, then, we are not really 
free to choose any gauge we wish. The choice is imposed by our observational 
situation and by the data we use in determining the best-fit FLRW model. 
Thus, though it may be may still be important to identify and use gauge- 
invariant quantities to resolve some issues, the perturbed quantities specified 
by the observational situation are observationally based and observable, and 
thus will reflect both the deviations of our space-time from FLRW and the 
errors and lacunae in the data used to determine our best-fit FLRW model. 
With regard to these latter, we are automatically led to modelling these as 
the imprecision or uncertainty in our identification of the background FLRW 
space-time using the data we have available. 
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5 Perturbed null Raychaudhuri equation 



Our first step in solving the perturbed spherically symmetric field equations 
is the solution of the perturbed version of the null Raychaudhuri equation 
( 25|) on our light cone w = wq to find the redshift z = z(y). All of our 
observational data on C~(po) is given as a function of z. In order to proceed 
with the integration we need to find z = z(y), so that we can give the 
data as a function of y for w = wq. Following Stoeger et al, (1992c), since 
C(wo,y) = t , the exact null Raychaudhuri equation fl25|) can be written as 



r 



r A 2 



and after some manipulation put in the form: 



d_ 

dy 



z— (l + z) 
dz 



1 dz 
--A r (l + z)M (z) — 

2 dy 



(37) 



(38) 



This equation has a first integral: 



dz A ( dr 0p dr 0+ 
dy V dz dz 



1 



-i 



'l + z) 



x M - - J (1 + z)(r 0F + r 0+ )(M 0F + M 0+ )dz 



(39) 



where we have written the data for the spherically symmetric perturbation 
problem in the form 



r Q (z) = r 0p (z) + r 0+ (z) 



(40) 



M (z) = M 0f (z) + M 0+ (z). (41) 
Here the subscripts F and + denote the FLRW and non-FLRW compo- 



nents of the data respectively. It is important to point out ( [[Stocgcr (1987) 



ttllis and Stoeger (1987)| , [Stoeger et al. (1992c)|| that in the case of FLRW 



both ro F and Mq f have very particular functional forms and if observations 
cannot be fit by those functional forms the universe is not FLRW. 

It was shown by Stoeger et al, (1992a) that considering r 0p and M 0f 
as zeroth-order quantities and r and M 0+ as first-order quantities, we can 
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write equations for the successive orders in the perturbation of dz/dy as 
follows. Zeroth-order equation: 



;i + ^)" 2 x ll-^j(l + z)r 0F M 0F (z)dz\ (42) 



First order equation: 




-A 



dr 0+ /dz 
(dr 0F /d2;) 2 



dr. 



-i 



Of 



dz 



[1 + zy 



(1 + z)- 2 <h - - / (l + z)r 0F {z)M 0F {z)dz 



x < -- J (1 + z)(r 0F (z)M 0+ (,2) + r 0+ (z)M QF (z))dz 



(43) 



Integrating equations ( f42|) and fl43|) yields 

= 2!f(|/) 21 = 21 (y), 

and to first-order we write 



(44) 



z iv) = z F{y) + z 1 {y). 



(45) 



Obviously, in the same way we can go to higher orders in the perturbation 
series. In this paper, we restrict ourselves to first order. Now that we have 
shown how to determine the redshift Z clS Bb function of y from the data 
functions, we can proceed to find the perturbed metric functions in terms of 
y and w. 

Before doing so, however, it is helpful to recognize that the perturbed 
data functions r 0+ (z) and M 0+ (z), as well as their FLRW background coun- 
terparts, will themselves be constructed from discrete data - with their many 
gaps and errors - by fitting some continuous function to them, a power law, 
or a polynomial, for instance (||Stoeger et al. (1992a)||). More precise or more 



complete data will necessitate a new fitting, to obtain improved data func- 
tions. The amount of uncertainty in the solutions due to gaps and errors 
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in the data will then be represented by a component of the perturbed met- 
ric functions. These can be tracked separately from the genuine deviations 
from FLRW simply by labelling the errors in the data functions by the usual 
indicators (e. g. error bars) and then determining the uncertainty or im- 
precision these errors induce in the solutions to the perturbation equations. 
The genuine deviations from FLRW will be those deviations which are larger 
than the error bars. These, too, can be traced in the same way through the 
perturbation equations to see whether or not they determine a model which 
is 'almost FLRW. 

Finally, it is important to recognize that there is a limit y* to the null co- 
moving radial coordinate out to which reliable data can be obtained [Stoeger 
et al. (1992a)]. The initial data on C~(po) for greater values of y than 
will remain unknown. Consequently, we can really only perform our in- 
tergrations for y < y*. Beyond that we really do not have the observa- 
tional data needed to determine the geometry of space-time. We do, how- 
ever, have CMWBR data from a redshift of z « 1500 ( jStoeger et al. (1995 



[Stoeger et al. (1997)|| ). However, that really does not serve to constrain ad- 
equately the space-time at more recent epochs and at less than those largest 
cosmological length scales. Note, however, that out to we can, in this 
spherically symmetric case, both predict and retrodict to all allowed values 
of w ( [ [Stoeger et al. (1992a)| | ) . This is because of the assumed spherical sym- 
metry, which effectively converts the usual hyperbolic equations of general 
relativity (for which prediction on the basis of given data is, strickly speak- 
ing, impossible) to ordinary differential equations for which these predictions 
to the future are legitimate. 



6 Completing the Integration Procedure 

In order to complete the integration of the perturbed field equations we 
introduce an improvement of the integration scheme recently developed by 
Araiijo and Stoeger (1999) for the exact spherically symmetric case that 
simplifies our present task. 

This modification of the integration scheme can be described schemati- 
cally as follows: 

(i) Having solved the null Raychaudhuri equation one finds A(wo, y) from 
equation ([T7|) . 

(ii) Since y is chosen to be a comoving radial coordinate the functional 
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dependence of A(w,y) with respect to y can not change as we move off our 
light cone. 

(iii) We have mentioned above the freedom of rescaling the time coor- 
dinate w that is effected by choosing A(w,0). So, given (i) and (ii) above, 
this freedom effectively corresponds to choosing the functional dependence 
of A(w, y) with respect to w in any way we like constrained only by the form 
of A(w ,y). In our expression for A(w ,y) is hidden an implicit dependence 
on w. We need to extract that dependence and make it explicit, so that we 
can then determine the general dependence of A on w and proceed with the 
integration. In general, this is not simply achieved by replacing wq with w 
because besides the w dependence arising from setting w = w when we 
write equation (|8|), there may be another part of the w dependence which 
derives from integration constants of the null Raychaudhuri equation and 
remains through the entire problem. With the aim of clarifying these com- 
ments we take for instance the expression for A(wo,y) that is obtained from 
(|T7| ) and the solution to the null Raychaudhuri equation with FLRW k = 
data ( fAraujo and Stoeger (1999)1 ), that is, 

Mw ,y) = 1 ^(i- 1 -) 2 (46) 

H w \ w J 

where Hq = H(wq, y) is the Hubble parameter. At this point assume that 
we arbitrarily set the w dependence for A and proceed with the integration. 
The next step is then the solution of equations ( |28| ) and fl30|) to determine B 
and C respectively. Their general solutions are: 

B = - J A'dw + l(y) (47) 
where l(y) is determined from the condition A(w ,y) = B(wo,y), and 



C 



^(-u, ) 1/2 / Achv + /,(//) 



2/3 

(48) 



where h(y) is determined from the data r = C(w , y). 

Now, we note that examining the central conditions (|15D on C(w, y) force 
B(w,y) to be of a certain form regarding its w dependence, which in turn, 
through equation (f4~7|) constrain A(w, y) to have a definite form. In the 
present example it is clear that unless 

4»,V) = ^~ (49) 
H w \ w J 
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the central conditions and the form A(w , y) given by the data and the solu- 
tion of the null Raychaydhuri equation will not be satisfied. 

(iv) B(w,y) and C(w,y) are then determined by integrating equations 
(|28| ) and (|30|) with respect to w. The arbitrary functions of y that arise from 
these integrations are determined by the conditions A(wo, y) = B(wq, y) and 
C(w ,y) = r (y) respectively B(w,y) and C(w,y) are further constrained 
by the fact that they have to satisfy the central conditions fllBT) . Now, it 
is clear from an examination of these equations that unless A(w,y) has a 
very specific functional dependence on w the resulting solutions B(w,y) and 
C(w,y) will not satisfy the central conditions. That implies that, although 
we can find solutions to the field equations, it does not guarantee that the 
null surface on which we assume we have the data is a past light cone of our 
world line ( ||Ellis et al. ( 1985^11 ) . So we conclude that (i), (ii) and the central 
conditions flip] ) remove the freedom of rescaling the time coordinate w and 
completely determine A(w,y). Thus, all the coordinate freedom in y and w 
has been used up, and this, together with the determination of the best-fit 
FLRW model to the data serves to specify the gauge (see Section 4). 

We now proceed to show how we can solve the perturbed spherially sym- 
metric field equations following the procedure described above. 

Combining equations (§) and ([|5|), we obtain to first order: 



l + z F (y) [l + z F (y)}< 



(50) 



Equating orders separately yields 



1 + z F (y) 

and 



R(w ,y)=> A = R (51) 



Z 2 (w , y) = 2A R(w , y) . (52) 

[l + ZF{y)\ 2 

In principle, the freedom of rescaling the time coordinate w together 
with our choice of a comoving radial coordinate y means that we can set 
the time dependency of Z 2 in any way we like constrained only by the form 
of Z 2 (wo,y). So at this stage, given this choice, we already have Z 2 (w,y), 
which is fundamental for proceeding with the integration of the remaining 
field equations for the determination of P 2 (w,y) and h^iw.y). But, as we 
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discussed above, this freedom does not really exist because of the further 
constraints imposed by the central conditions. So in order to proceed with 
the integration we have to set the correct time dependency of Z 2 at this 
point. 

Substituting ( fL8|) and ( |19|) into fl28f ) and collecting the first order terms 
yields 



[r) 2\r) 2{R / 
Integrating ( |53| ) with respect to w we find 



(53) 



(3 2 1 f fZ 2 \' 1Z 2 . . 
R = 2J [-I) dw -2-R + 3{y), 



(54) 



where g(y) is determined by the condition (3 2 (w ,y) = —Z 2 (w ,y). 

In order to complete the solution, we have to find h,2 2 (w,y) from the 
perturbed form of equation (|30|). But first we must calculate the "potential" 
W(y) and u (y) from equations (|29D and (|32|) , respectively. That in turns 
involves the calculation of h 22 (wo, y). To do this we must solve the perturbed 
form of equation ( [27D on our past light cone w = w . Using equations fllT 
( ]TP| ) and (PU|), the first order perturbation equation is given by 



\Rf 

h 22 



+ 



(Rfy (h 22 \ (Rfy fh 22 



2R 3 f 



2Rf \Rf 

(Rjy (Rjy + r 



2Rf \Rf, 







2R 





R 2 h 22 _ (Rjy /U' h 22 [(Rf) 
4R 3 P 2Rf \Rf) 4i? 3 / 3 



+ 



1 f(3 2 



(R}y + + gKgffi? + ^tW)'] 2 



R \ R J 2Rf 2R 3 f 2R 3 f 

Since on our light cone /3 2 (w , y) = — Z 2 (u7 , y) we obtain at w = w 



(55) 



'^22' 



+ 



(Rfy (k 22 

Rf \Rf, 



(Rfy A22V 
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'22 



R 3 f 3 
R [ R 



(Rf)(Rf)' + R 

(RfY - 







z 2 




R3 ( 


22 


(RfY ( 


'h 22 



\/l2 



Rf \Rf, 



2R 3 f 3 



+ 2 -n\RfY-^. 

R\RJ y ' Rf 



(56) 



This is a standard linear equation whose integrating factor can be easily 
calculated and yields the solution 



' h 2 2 

.Rf, 



{w ,y) 







RfJ 


r{- 


+R 


/ ^22 

\Rf 


R 2 h 22 


2R 2 f 2 




'9" 



+ 



'22 



R 2 P 



(Rf) 



( h<12 
\Rf, 

a* 



, f h 22 \ + h 22 [(Rf)f 



(RfY (RfY 



Rf, 



2R 2 P 



(Rf)'f - Z 2 dy. 



(57) 



We can now use these results to construct the potential W(y) given in 
equation (^). Note that knowing C', C, B and A at any one time, for 
instance w = Wq, determines W(y) for all times since it is independent of w. 
In our case here we have just determined h 22 (wo, y). From the data we know 
h' 22 (wo,y) and the from the solution to the null Raychaudhuri equation we 
have Z 2 (wo,y) = —/3 2 (wo,y). We can then write (evaluating at w = Wq) 



W(y) = W F (y) + W 1 (y) 



(58) 



where, 



W F (y) 
W x (y) 



1 

R 

Z 2 



2R 3 



[(RfY + (RfY] 
[(RfY - (RfY] 



B 2 - 1 



h 22 



+ 



h 22 

Rf, 



(59) 
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At this stage, since we have already solved the null Raychaudhuri equation 
in Section 5 to find z = z(y), we can write u Q given in (p4]) as a function of 
y to first order as 

u (y) = u F (y) +wi(y) (60) 

where, 



W , 



w 0l 



+ ±-(Rf)'(Rf)W F (l + z F ) - ^{Rf)[{Rf)'?{l + z F f 



h22 1 n , \ 
4R/ 



{Rf)\Rf)W l + \(^ (Rf)W F 



+ ±-(Rf)' (Rf)W FZl ' " - - 



2AI 



;i + z F f 



x (i2/)[(i2/)'] 2 ' 



(61) 



Substituting (p0[), (|58|) and flBOP into ( |3~0|) and equating orders separately 
we obtain the following first order equation for /i 2 2 



" W^Rl = R ( 2WpWl ~ 2 "° l/i?/ ) " ^ (62) 

where F = (W| - 1 - 2uj 0F /Rf) . Equation (H) is a standard linear 
equation whose general solution is 



1 



e a(y)w 
1 

p a(y)w 



e a{y)w b(y)dw + c{y) 

a(y)w 



Kv) 



+ c(y) 



a(y) 



+ e' 



a(j/)ru\-l 



(63) 
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where, 



a(y) 



U0 F 



F-\y) 



RP 



b(y) = R(2W F W l -2uj Q jRf)F-\y)-^-F{y) 



(64) 



and c{y) is a function determined by the data /^("^cb y)- 

To summarize, once we solve the null Raychaudhuri equation we find 
Z 2 (wq, y) from equation (0), then according to our discussion in the begining 
of this section concerning the use of the central conditions we determine 
Z 2 (w,y). We then proceed to find P(w,y) (equation (|54T)). Finally h 2 2(w,y) 
(equation (jBHj) ) is determined completing the solution. 

This procedure gives us the full solution to the first-order perturbed equa- 
tions in terms of data functions for number counts and observer area distance 
on our past light cone. As we have already mentioned briefly in Section 5, we 
could go on in a similar way to solve the perturbation equations for second 
order and higher. In order to determine whether or not our perturbation 
treatment is valid, we shall have to look at higher-order solutions. Does the 
perturbation series we obtain converge, or at least give some evidence of con- 
verging, for the data functions r 0+ (z) and M 0+ (z) we input? If it does, and if 
those solutions remain convergent for all values of w, then we can affirm in a 
rigorous sense that our universe is 'almost FLRW on the scales represented 
by the data. If they are not, then the perturbation treatment is invalid, and 
our universe is not 'almost FLRW on those scales. 

Of course, in doing this, we must also track separately the 'error com- 
ponent' and the 'genuine perturbation' component of the metric deviations 
from FLRW. We briefly indicated at the end of Section 5 how this can be 
done. 

In the next paper, we shall present the solution to the general FLRW 
perturbation problem for observational data on our past light cone. That 
case will be similar in many ways, but much more complicated in others. 
In fact, it is only with the insights gained from this detailed study of the 
solution to the spherically symmetric case that we see how to proceed with 
the more general case. 
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